Substellar Companions to Main Sequence Stars: No Brown Dwarf 

Desert at Wide Separations 
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ABSTRACT 

We use three field L and T dwarfs which were discovered to be wide companions to known stars 
by the Two Micron All-Sky Survey (2MASS) to derive a preliminary brown dwarf companion 
frequency. Observed L and T dwarfs indicate that brown dwarfs are not unusually rare as wide 
(A > 1000 A.U.) systems to F-MO main-sequence stars (M > 0.5Mq, My < 9-5), even though 
they are rare at close separation (A < 3 A.U.), the "brown dwarf desert." Stellar companions 
in these separation ranges are equally frequent, but brown dwarfs are > 10 times as frequent for 
wide than close separations. A brown dwarf wide-companion frequency as low as the 0.5% seen 
in the brown dwarf desert is ruled out by currently-available observations. 

Subject headings: binaries: general — stars: low-mass, brown dwarfs 



1. Introduction 

Understanding the processes, and distinctions 
between, star formation, binary formation and 
planetary formation is a major goal of stellar as- 
tronomy. A necessary step towards that goal is to 
understand the frequency of "brown dwarf" and 
"planetary" companions as a function of separa- 
tion, primary mass, and secondary mass. It is 
now well-known that brown dwarfs are very rare 
as close companions to F-M dwarfs (the "brown 
dwarf desert," Marcy & Butler 2000) yet they are 
quite common in the field (Reid et al. 1999) and 
in open clusters (Bouvier et al. 1998; Martin et al. 
1998). 
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The purpose of this letter is to show that cur- 
rent data allow an estimate of the wide (A > 
1000 A.U.) brown dwarf companion frequency to 
near-solar mass main sequence stars, despite the 
present lack of well-defined searches for wide com- 
panions. The observational constraints and re- 
sulting wide companion frequency are discussed in 
Section 2 and the differences between this fraction 
and the close companion frequency are discussed 
in Section 3. 

2. Observational Constraints 

2MASS searches for ultracool dwarfs are based 
on colors and magnitudes; although designed to 
find isolated field dwarfs, they are not biased 
against widely separated (> 40 arcsec) compan- 
ions whose photometry is uncontaminated by the 
primary star. Three published brown dwarf com- 
panions with separations A > 1000 A.U. have 
been identified in the course of 2MASS searches 
for isolated field brown dwarfs. The L4.5 dwarf Gl 
417B has an estimated mass of 0.035 ± O.O15M0 
and age 0.08-0.3 Gyr while the L8 dwarf Gl 584C 
has a mass of 0.060 ± 0.015Mo and an age of 1.0- 
2.5 Gyr (Kirkpatrick et al. 2000, 2001). The T 
dwarf Gl 570D has a mass of 0.050 ±0.020Mq and 
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an age of 2-10 Gyr (Burgasser et al. 2000). All 
are definitely substellar via spectroscopic crieria 
as well, since both L dwarfs pass the lithium test 
while the T dwarf is too cool to be a star. All 
three have been confirmed as companions by their 
common proper motion. These reflect the results 
from 2MASS searches of only a fraction of the sky 
— wide brown dwarf companions continued to be 
discovered. Wilson et al. (2001) will describe the 
discovery of three additional L dwarf secondaries 
at 880 A.U., 1090 A.U., and 2460 A.U. 

The usual procedure for determining compan- 
ion frequency is to survey a number (N) of poten- 
tial primary stars. If n brown dwarfs with spec- 
tral type L are found, the frequency of L dwarf 
companions (in the searched separation range) is 
simply /l = In the case of imaging surveys, 
only a fraction (j/^) of the brown dwarfs can be de- 
tected as L dwarfs, since brown dwarfs fade to very 
low temperatures and luminosities. Modelling of 
Ul allows an estimate of the true brown dwarf fre- 
quency to be fhd = . In principle, the Two Mi- 
cron All-Sky Survey (2MASS) allows a search for 
wide companions to nearby stars, allowing and 
fbd to be determined. A complete search, however, 
has not yet been made, due to the complexities of 
the task. 

It is nevertheless possible to estimate the fre- 
quency of wide-separation brown dwarf compan- 
ions, since the frequency can also be expressed as 

/. Pcomp Pcomp Pbd 9wPbd 9wPL 

Jbd = = = = 

Pstar Pbd Pstar Pstar VLPstar 

In this equation, pstar is the space density of 
stars (the potential primaries), Pcomp is the space 
density of wide brown dwarf companions, Pm is 
the space density of field brown dwarfs found by 
2MASS, and is the frequency of field brown 
dwarfs that have a stellar primary at A > 1000 
A.U. Each of these variables can be estimated from 
published data. Separations of 1000 A.U. corre- 
spond to 40 arcseconds at 25 parsecs distance, 
and larger values at closer distances. For these 
wide systems, the presence of the main sequence 
star does not affect the identification of the brown 
dwarf in searches for isolated field objects. In the 
final step, we have assumed that for all brown 
dwarfs is the same as for L dwarfs. ^ 

Some of the field L dwarfs used in our calculation will 



The observed space density of L dwarfs can 
be determined using the Kirkpatrick et al. (1999) 
sample of L dwarfs. Adopting the limiting magni- 
tude Kg = 14.7 and the Kirkpatrick et al. (2000) 
spectral type-Mi^^ relation, the Schmidt (1968) 

— technique yields an L dwarf space density 
oFpL = 0.0057 ± 0.0025 pc"'\ The space den- 
sity of T dwarfs is highly uncertain but compa- 
rable (Burgasser et al. 1999). Only a small frac- 
tion (ul) of brown dwarfs are observable as L, or 
even T, dwarfs. We can estimate j/i in two ways. 
// the mass function of companions is similar to 
that of isolated field brown dwarfs, then we can 
use the correction determined for isolated brown 
dwarfs. Reid et al. (1999) estimate that if the sub- 
stellar mass function can be described as a power 
law ^ oc extending down to 0.01 A/q, then 

a w 1.3 and pbd ~ 0.10 b.d. pc~^ with large un- 
certainties. For the cases a = 0.0, 0.6 and 1.0, 
then Pbd « 0.02, 0.04 and 0.07 b.d. pc~^ respec- 
tively. Surveys of the Pleiades are consistent with 
a « 0.6 - 1.0 (Bouvier et al. 1998; Martin et al. 
1998), and neither that cluster nor the field is con- 
sistent with a = 0. We adopt ptd = 0.07 ± 0.03 
b.d. pc~^; hence, the ratio of L dwarfs to brown 
dwarfs is w 0.08. 

Alternatively, we may attempt to estimate the 
parameter directly from the properties of the 
L dwarf companions. Gl 584C is just at the lim- 
its of detectability; since its age is in the range 
1.0-2.5 Gyr, it is only detectable for 0.1-0.25 of 
the age of the disk. Gl 41 7B is easily detectable 
as an L4.5, but would have been detectable down 
to spectral type L8. Comparison to evolutionary 
tracks indicates it would then be visible to an age 
of 0.3-1.0 Gyr. The companions themselves then 
imply a correction factor in the range 0.07 - 
0.18 — consistent with the isolated brown dwarf 
estimate. 

The fraction of the "isolated" field brown 
dwarfs which are actually wide companions to 
main sequence stars can be estimated from the 
2MASS field brown dwarf searches. Over 100 
field L dwarfs are now published, primarily from 
the 2MASS (Kirkpatrick et al. 1999, 2000), DE- 
NIS (Delfosse et al. 1999), and SDSS (Fan et al. 

not be brown dwarfs; however, since both Gl 417B and Gl 
5840 are confirmed brown dwarfs, accounting for this ef- 
fect will only increase the derived brown dwarf companion 
frequency. 
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2000) surveys, of which two are wide compan- 
ions to near-solar mass stars. Many of the L 
dwarfs, however, are at large distances (> 25 pc), 
where it would be difficult to discover an L dwarf 
companion due to proximity of the primary star. 
We choose to consider only L dwarfs within 25 
parsccs because at this distance, we can be con- 
fident that bright (F-K) primaries are cataloged 
and the L dwarf companions at A > 1000 A.U. 
will be detectable. In this case, 2 of 40 L dwarfs 
(Kirkpatrick et al. 2000, Table 4) have primaries, 
indicating the fraction of field L dwarfs which 
are actually wide companions is = 0.05 ± 0.04. 
This estimate is supported by the ongoing 2MASS 
T dwarf search (Burgasser et al. 1999), which has 
now discovered 17 field T dwarfs Burgasser (2000). 
Since T dwarfs are simply the older counterparts 
of L dwarfs, the fraction of 'field' dwarfs which 
are companions should be roughly the same. One 
(Gl 570D) of the T dwarfs is the fourth member of 
a K/M/M/T quadruple system (Burgasser et al. 
2000), leading to an estimate of = 0.06 ± 0.06. 
(Note that Gl 229B is too close to Gl 229A to 
be counted or even detected by 2MASS.) Com- 
bining the L dwarf and T dwarf values, we adopt 
gyj = 0.05 ±0.03. 

The space density of stars which are poten- 
tial primaries for the brown dwarfs can be deter- 
mined from the nearby star catalog as modified 
by Hipparcos parallaxes. All of the primaries to 
wide L/T dwarf companions (including the new 
Wilson et al. 2001 systems), have My < 9.5 (or 
Mass > O.5M0). Jahreifi & Wielen (1997) find 
that the space density of stars with My < 9.5 is 
0.020 ± 0.001 pc-3. These are F to MO main se- 
quence stars. We do not consider fainter primaries 
becaiisc the incompleteness of the late-M dwarf 
catalog means that we might not recognize their 
presence near an L dwarf. We therefore have in- 
sufficient information to calculate the M/L dwarf 
binary statistics. * We can now estimate the 
wide brown dwarf companion frequency to main 
sequence stars. Taking the stars with My < 9.5 as 

'Note that the early-M dwarfs with 9.5 < My < 13.5 con- 
tribute another 0.039 ± 0.008 stars pc~'^, implying that we 
could decrease the derived brown dwarf frequency by a fac- 
tor 2.9 — but in this case we face the puzzling fact that the 
chance of all 5 primaries having My < 9.5 is only ~ 0.4%. 
An analogous difficulty in interpretation would occur if a 
search for companions amongst F-M dwarfs found compan- 
ions only around the F, G and K dwarfs. 



the available primaries, we find /l = 0.014±0.011. 
These, however, represent only a small fraction of 
the total brown dwarf population. We estimate 
fbd = 0.18 ± 0.14 by assuming a mass function 
a = 0.7; the yL values from the companions them- 
selves suggest values in the range 0.08-0.20 ±0.14. 

The large uncertainties — due to the fact that 
the estimates are based only 3 companions — in- 
dicate that a larger survey is needed to determine 
the brown dwarf companion fraction; we certainly 
cannot prefer fbd = 0.18 over fbd = 0.05 given 
the uncertainties. However, it is important to re- 
alize that the error bars are non-Gaussian and 
that very low companion fractions (fbd) are al- 
ready ruled out with statistical significance. The 
situation is illustrated in Figure 1. We run Monte 
Carlo simulations of samples of 57 brown dwarfs 
and determine the percentage that have at least 3 
wide companions. (For each value of fbd, we run 
200,000 simulations.) The solid line plots our pre- 
ferred scenario: A brown dwarf space density of 
0.07 pc~^ (equivalently, i/l = 0.08) and all pri- 
maries with My < 9.5 considered. It is possi- 
ble to consider other scenarios. If pbd = 0.02, 
equivalent to assuming a = or = 0.29 (long- 
dashed line), then values of fbd as low as 0.015 
can still be excluded at the 95% confidence level. 
This corresponds to a scenario in which nearly all 
brown dwarfs are relatively massive and therefore 
detectable as L dwarfs for many Gyr; in this case, 
the observed L dwarfs represent a larger fraction 
of the total brown dwarf population. We finally 
consider an extreme model, in which pbd = 0.02 
and Pstar = 0.059 pc^"*. In this case, the brown 
dwarf binary fraction can be much lower; in other 
words, we can reduce the reported brown dwarf 
fraction by averaging the many brown dwarf com- 
panions to M > 0.5M© primaries with the more 
numerous population of lower mass stars for which 
no wide brown dwarf companions have ever been 
detected. 

3. Discussion 

The existence of the "brown dwarf desert" sep- 
arating stars from "planets" at separations less 
than 3 A.U. is well established (Halbwachs et al. 
2000). The review of Marcy & Butler (2000) finds 
that less than 0.5% of F-M dwarfs have brown 
dwarf companions down to 0. 01 Mq on the basis of 
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over 500 stars. The Mazeh et al. (1992) analysis of 
the 164 nearby G dwarfs studied by Duquennoy & 
Mayor (1991) indicates that 13 ± 3% of G dwarfs 
have stellar companions in this separation range, 
and that the mass distribution of these close com- 
panions is lacking in very-low-mass stars (a "red 
dwarf steppe" ?) compared to both wider compan- 
ions and the field mass function. On the basis of 
Duquennoy & Mayor (1991)'s analysis (their Fig- 
ure 7), we estimate that 12±3 % of G dwarfs have 
stellar companions with separation A > 1000 A.U. 

Our evidence that F-MO dwarfs have a signifi- 
cant {fbd = 18 ±14%) population of very- low- mass 
companions at wide separations indicates that 
the orbital separation distribution of brown dwarf 
companions is not simply a scaled-down version of 
the stellar companion distribution. G dwarfs have 
approximately equal numbers of stellar-mass com- 
panions at separations A < 3 A.U. and A > 1000 
A.U. (riooo/3 ~ 0.9); in contrast, we estimate 
their brown dwarf companions are at least four, 
and probably many more, times more common 
at the larger separations (r"iooo/3 ^ 4). Within 
the very large uncertainties, the number ratio of 
stellar and brown dwarf companions may be near 
unity, similar to the ratio for isolated stars and 
brown dwarfs (Reid et al. 1999), although a con- 
siderably larger sample is needed to investigate 
this. The fact that all the primaries are rela- 
tively massive suggests that wide brown dwarf 
companions do not form around — or are not 
retained by — less massive (M < 0.5Mq) pri- 
maries, although the incompleteness of the nearby 
star catalog for M dwarfs may contribute to this 
effect. The extreme limit to this is noted by Reid 
et al. (2001), who find that L dwarf primaries 
lack companions beyond ^ 10 A.U. Where the 
"brown dwarf desert" for F-MO dwarf primaries 
ends is unclear at present, but searches around 
stars in the range 1-100 A.U. have generally had 
little success even though this range is the peak of 
stellar companion distribution for both G and M 
dwarf primaries (Duquennoy & Mayor 1991; Fis- 
cher & Marcy 1992). Most recently, Schroeder et 
al. (2000) failed to discover any brown dwarfs de- 
spite a sensitive Hubble Space Telescope search at 
separations of 1 — 60 A.U. around 23 stars, while 
the extensive, highly sensitive search of 107 stars 
by Oppenheimer et al. (2001) found only 1 brown 
dwarf in range 40 - 120 A.U. (Note that these 



voliurie- limited samples are dominated by M dwarf 
primaries with M < 0.5M©.) The fraction in the 
range 100-1000 A.U. cannot yet be constrained, 
but the recent discoveries of three L brown dwarf 
companions at those separations (Rebolo et al. 
1998; Goldman et al. 1999; Gizis et al. 2001) sug- 
gest that they may be as common as wide compan- 
ions and that 2MASS should detect many more. 

While the uncertainty of our derived compan- 
ion fraction is large, perspective on the brown 
dwarf desert may be gained by considering the 
numbers of stars and brown dwarfs within 25 par- 
sees. There are 1297 main sequence stars with 
My < 9.5 in this volume according to the Jahreifi 
& Wielen (1997) luminosity function. Radial ve- 
locity surveys indicate that they will have fewer 
than 6.5 brown dwarf companions within 3 A.U. 
In constrast, three wide-separation (two L, one 
T) brown dwarf companions are already known 
to such stars even though 1) only a small frac- 
tion of the sky has been searched, 2) T dwarfs are 
not detectable all the way to 25 parsecs, and 3) 
older/lower-mass brown dwarfs are undetectable. 
Our numbers suggest another ^ 16 ± 9 wide L 
dwarf companions within 25 parsecs will be found, 
representing only the young, massive tip of the 
brown dwarf population. 

Combined with the Reid et al. (2001) observa- 
tion that 20% of L dwarfs have H.S.T. resolved 
(> 0.1 arcsec) brown dwarf companions within 10 
A.U. but that wider companions are lacking, (see 
also Koerner et al. 1999), there is strong evidence 
that the frequency of brown dwarf companions is 
strongly dependent on both primary mass and or- 
bital separation. There is a "brown dwarf desert" 
at < 3 (or < 100?) A.U. for F to mid-M main 
sequence primaries, and another desert for wide 
(> 20 A.U.) brown dwarf doubles. This situation 
arises naturally in some theories of star formation 
(Boss 1986; Bate 2000). 

4. Summary 

We use the field L and T dwarfs which were 
discovered to be wide companions to known stars 
to derive a preliminary brown dwarf companion 
frequency. The observed L and T dwarfs indicate 
that brown dwarfs are not unusually rare as wide 
(> 1000 A.U.) systems, even though they are rare 
(< 0.5%) at close (< 3 A.U.)separations. The 
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current data indicate that ^ 1% of My < 9.5 pri- 
maries have a wide L dwarf companion; the brown 
dwarf fraction should be substantially (5—13 
times) higher. Stellar companions in these sep- 
aration ranges are equally frequent. 

Our estimates indicate that continued searches 
for wide brown dwarf companions should not be 
discouraged by the existence of the "brown dwarf 
desert" at close separations. There is now strong 
evidence — from both the wide companions ex- 
amined in this work, and the close double brown 
dwarfs examined by Reid et al. (2001) — that 
the brown dwarf companion frequency is a strong 
function of both primary mass and separation. A 
search using 2MASS for both stellar and brown 
dwarf companions with separation A > 100 A.U. 
will be quite rewarding. Furthermore, provided 
the older and lower-mass wide brown dwarf com- 
panions arc not stripped in the Galactic disk, the 
majority of the brown dwarfs should be intrinsi- 
cally cooler than Gl 570D, suggesting that a deep 
wide-field SIRTF search near G and K dwarfs is 
a promising avenue to extend our knowledge of 
brown dwarfs. 
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Fig. 1. — The percentage of samples in which 
3 wide companions are found in a sample of 57 
field brown dwarfs. The solid curve is our pre- 
ferred model {phd — 0.07,psiar = 0.020); the long- 
dashcd line is a model with few brown dwarfs 
(Pbd = 0.02,pstor = 0.020); the short-dashed curve 
is a model with few brown dwarfs and M dwarf 
primaries (pM = ^S)2,pstar = 0.059). The posi- 
tion of the brown dwarf desert (/bd < 0.005) and 
the 95% confidence level is also marked. 
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